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ABSTRACT 

Two decades ago "transitional disks" described spectral energy distributions (SEDs) of T Tauri stars 
with small near-IR excesses, but significant mid- and far-IR excesses. Many inferred this indicated 
dust-free holes in disks, possibly cleared by planets. Recently, this term has been applied disparately 
to objects whose Spitzer SEDs diverge from the expectations for a typical full disk. Here we use 
irradiated accretion disk models to fit the SEDs of 15 such disks in NGC 2068 and IC 348. One 
group has a "dip" in infrared emission while the others' continuum emission decreases steadily at all 
wavelengths. We find that the former have an inner disk hole or gap at intermediate radii in the 
disk and we call these objects "transitional" and "pre-transitional" disks, respectively. For the latter 
group, we can fit these SEDs with full disk models and find that millimeter data are necessary to 
break the degeneracy between dust settling and disk mass. We suggest the term "transitional" only 
be applied to objects that display evidence for a radical change in the disk's radial structure. Using 
this definition, we find that transitional and pre-transitional disks tend to have lower mass accretion 
rates than full disks and that transitional disks have lower accretion rates than pre-transitional disks. 
These reduced accretion rates onto the star could be linked to forming planets. Future observations 
of transitional and pre-transitional disks will allow us to better quantify the signatures of planet 
formation in young disks. 

Subject headings: accretion disks, stars: circumstellar matter, planetary systems: protoplanetary disks, 
stars: formation, stars: pre-main sequence 



1. INTRODUCTION 

Disks around T Tauri stars (TTS) are thought to be 
the sites of planet formation. However, many questions 
exist concerning how the gas and dust in the disk evolve 
into a planetary system and observations of TTS may 
provide clues. There are some objects in particular that 
have gained incr easing attent i on in this regard. Over 
two decades ago, IStrom et~a l. (1989) detected "possible 
evidence of changes in disk structure with time" as evi- 
denced by "small near-IR excesses, but significant mid- 
and far-IR excesses" indicating "inner holes" in disks. 
Those authors proposed that these objects were 11 in tran- 
sition from massive, optically thick structures that ex- 
tend inward to the stellar surface, to low-mass, tenuous, 
perhaps post-planct-building structures." 
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Usage of the term "transitional disk" gained sub- 
sta ntial moment u m in th e literature after it was used 
by ICalvet et al.l ((2005bf ) to describe disks with in- 
ner hol es using data fro m the Spitzer Space Tele- 
scope's (|Werner eF al. 2004) Infrared Spectrograph (IRS; 
iHouck et al.l |200*4TT Before Spitzer, the spectral en- 
ergy distributions (SEDs) used to identify disks with 
inner holes were based solely on near-infrared (NIR) 
ground-based photometry and IRAS mid-IR (MIR ) 
photometry ()Strom et all 119891: ISkrutskie et~aTl 119901) . 
Spitzer allowed the opportunity to study these objects 
in greater detail. The unprecedented resolution and 
simultaneous wavelength coverage (~5 and 38 /iim) of 
Spitzer IRS unco vered new details r egarding these disks 
(ID'Alessio et al.l 120051: ICalvet et ai]|2005bt iFurlan et all 
|2006J). Some SEDs had nearly photospheric NIR (1- 
5 /im) and MIR (5-20 /im) emission, coupled with sub- 
stantial emission above the stellar photosphere at wave- 
lengths beyond ^20 /im. Others had significant NIR 
excesses relative to their stellar photospheres, but still 
exhibited MIR dips and substantial excesses beyond 
~20 /im. 

Detailed modeling of many of the above-mentioned 
SEDs has been performed. SEDs of disks with little or 
no NIR and MIR emission have been fit with models 
of inwardly truncated optically thick disks (jCalvet et al.l 
12001 120051* rEspaillat etafll2007bt l2008bh . The inner 
edge or "wall" of the outer disk is frontally illuminated 
by the star, dominating most of the emission seen in the 
IRS spectrum. In this paper, we refer to these objects 
with holes in their dust distribution as transitional disks 
(TD). Some o f the holes in TP are relatively dust-free 
(e.g. DM Tau: ICalvet et al.ll2005bl : lEspaillat et alj|2010h 
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while SED model fitting indicates that others with strong 
10/im silicate emission have a small amount of optically 
thin dust in their disk holes to explain this feature (e.g. 
GM Aun lCalvet et alJl2005rJ : lEspaillat et alJl2010f >. For 
SEDs with substantial NIR emission accompanied by a 
MIR dip, we can fit the observed SED with an optically 
thick inner disk separated by an optically thin gap from 
an optically thick outer disk (Espaillat et al. 2007a). 
Here we call these pre-transitional disks (PTD). Like the 
TD, we see evidence for relatively du st- free gaps (e.g. UX 
Tau A; lEspaillat et al.ll2007aL 120101 ) as well as gaps with 
some small, optically thin dust to expla in strong 10/zm 
silicat e emis sion features (e.g. LkCa 15; Espaillat et al. 
I2007al I2010D . For many TDs and PTDs, the trunca- 
tion of the outer disk has been confirmed with sub- 
millimeter and millimeter interferometric imaging (e.g. 
DM Tau, GM Aur, UX Tau A LkCa 15:IHughes et al.l 
[20071 l2009t I Andrews et all [20091 l20Tol [20l¥ as well as 
NIR imaging (i.e. LkCa 15; iThalmann et al.ll20Tol ). In 
a few cases, the optically thick inner disk of PTD has 
bee n confirmed using the "yeilin^p " of near-infrared spec- 
tra (|Espaillat et al] l2008al 12010) and near-infrared in- 
terferometry has confirmed that the i nner disk is small 
(|Pott et al.ll2010t IQlofsson et al.l[20Tl . 

The distinct SEDs of TD and PTD most likely signify 
that these objects are being caught in an important phase 
in disk evolution. Many researchers have posited that 
these disks are forming planets on the basis that cleared 
dust re gions are predicted by planet formation mod- 
els le.g.lPaardekooper fc Mellema| [200l IZhu et al.H 2011t 
iDodson- Robinson fc Salvkl 120111) . Recently, a potential 
protoplanet has be en reported in the p re-transitional disk 
around LkCa 15 dKraus et al.l 12011) as well as around 
T Cha ([Huelamo et al.l 120 31). Stellar companions ca: i 
also clear the inner disk (|Artvmowicz fc Lubowl I1994T ) 
but many stars harbori ng transitional disks are single 
stars (Kr aus et al.ll20lH) . Even in cases where stellar- 
mass companions have not been ruled out, the large holes 
and gaps observed are most likely evidence of dynamical 
clearing. Photoevaporation cannot explai n disks with 
large cavities and high mass accretion rates (jOwen et 
120 111 ) and dust evolution alone can not explain the sharp 
decreases in surface density seen in the SED and inter- 
ferometric visibilities. 

Given the potential link between disks with gaps and 
holes and planet formation, interest in transitional disks 
has grown. Some studies have focused on further un- 
derstanding the behavior of the currently known mem- 
bers in this class of objects and have discovered IR vari- 
ability ([Muzerolle et al.ll2010t lEspaillat et aT]|2011l Ell). 
Other studies have taken a broader approach, working 
towards expandi ng the known number of disks undergo- 
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2010; Luhman et al. 


2010; Merin et al.l 


2010: Curric & Sicilia-Aguilar 2011). 



As the literature on Spitzer observations of TTS ex- 
pands, so does the terminology applied to disks around 



9 "Veiling" occurs when an excess continuum (Hartigan ct al. 
1989) "fills in" absorption lines, causing them to appear signifi- 
cantly weaker than the spectrum o f a standard star of the same 
spectral type (Hartigan et al. 1989). The veiling observed in pre- 
transitional disks is similar to that observed in full disks where the 
veiling has been explained by em ission from the inner d isk edge or 
"wall" of an optically thick disk (Muzcrollc ct al. 2003). 



TTS. These disks are referred to as primordial, full, 
transitional, pre-transitional, kink, cold, anemic, ho- 
mologously depleted, classical transitional, weak ex- 
cess, warm excess, and evolved. However, these 
terms are not applied consistently. The issue is high- 
lighted with the discrepancies in the reported frac- 
tions of transitional disks and disk clearing timescales 
(iMerin et al.ll2010t ILuhman et alj|2010t IMuzerolle et al " 
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20101 ). In many cases, the data is the same; the dif- 
ferences arise from nomenclature. 

If the goal is to better understand disk evolution, it 
is important to look past the nomenclature and deci- 
pher the underlying disk structure we can infer from the 
observations, while paying special attention to the lim- 
itations of both the observations and the tools we use 
to interpret them, namely disk models. To this end, 
we chose a sample for our study which incorporated ob- 
jects whose SEDs are similar to those that have been 
referred to as "transitional" in the literature. We fo- 
cus on 15 disks in NGC 2068 and IC 348 with Spitzer 
IRS spectra. NGC 2068 an d IC 348 are older (~2 Myr 
and ~3 Myr, respe ctively; iFlahertv fc Muzerolld 120081 : 
ILuhman et al.ll2003l ) and more clustered star-forming re- 
gions with higher extinction than Taurus, where most 
detailed studies of transitional disks have focused. Ac- 
cording to our definitions stated above, our sample con- 
tains five TD, five PTD, and five objects with decreas- 
ing emission at all IR wavelengths (i.e. negative MIR 
slopes) which we classify as full disks. Out of the ob- 
jects reported by IMuzerolle et ah! ((20100 . this includes 5 
of the 15 transitional and pre-transitional disks in IC 348 
and 4 of the 6 transitional and pre-transitional disks in 
NGC 2068. One of the objects classified as a full disk in 
that work is classified as a PTD here. 

In Section 2, we present optical photometry and spec- 
troscopy, infrared spectra, and millimeter flux densities 
that we used to compile SEDs for our objects, derive 
stellar properties, and measure mass accretion rates. In 
Section 3, we fit the SEDs of our objects with irradiated 
accretion disk models and in Section 4 we discuss the 
limitations of the models and the observations, disk se- 
mantics, and the mass accretion rates of transitional and 
pre-transitional disks. 

2. OBSERVATIONS & DATA REDUCTION 

Before conducting modeling, we first compiled SEDs 
and mass accretion rates for the targets in our sample 
(Table [T]). To supplement the existing information in the 
literature, we collected optical photometry for all of our 
IC 348 objects, optical spectra for three of the IC 348 
targets and all of our NGC 2068 targets, and millimeter 
flux densities for two of our IC 348 targets. We also 
present infrared spectra for all the targets in this paper. 
Below we discuss the details of the observations and data 
reduction. 

2.1. Optical Photometry 

During November 19-23, 2011 we used the 4K CCD 
imager on the 1.3-m McGraw-Hill telescope^ of the 
MDM Observatory to obtain UVRI photometry of the 
IC348 stellar cluster. The FOV of the instrument is 

10 http://www.astronomy.ohio-state.edu/~ jdeast/4k/ 
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21'x21' (0.315" per unbinned pixel). For these obser- 
vations, we used 2x2 binning (0.62" per pixel). The 
flat-field, overscan, and astrometric calibration were per- 
formed using an IDL program written by Jason East- 
mar0 specifically designed for the 4K CCD imager. 
Since large electronic structures are not stable enough 
to be reliably subtracted, we did not apply corrections 
using two dimensional biases. The photometric calibra- 
tion of all images was carried out using the standard pro- 
cedure and the daophot and phot cal packages in IRAF, 
with standard stars selected from iLandoltl (|1992j ). Pho- 
tometry is presented in Table [2] 

2.2. Optical Spectroscopy 

In order to obtain mass accretion rate estimates (Sec- 
tion 13.1.11 and Figur e ED), we collected M IKE double- 
echelle spectrograph ([Bernstein et al.ll2003D data on the 
6.5m Magellan Clay telescope for all of our targets in 
NGC 2068 (FM 177, FM 281, FM 515, FM 581, and 
FM 618) and some of our objects in IC 348 (LRLL 21, 
LRLL 67, and LRLL 72). The observations for the 
NGC 2068 and IC 348 objects were taken on February 
10-11, 2007 and January 19, 2009, respectively. We used 
a slit size of 0.7" x 5.0" and 2x2 pixel on-chip binning 
with exposure times of 800-1200s. Data were reduced 
using the MIKE data reduction pipeline o 

2.3. Infrared Spectroscopy 

Here we present Spitzer IRS spectra for each of our tar- 
gets. Spectra for NGC 2068 and IC 348 were obtained 
in Program 58 (PI: Rieke) and Program 2 (PI: Houck), 
respectively. All of the observations were performed 
in staring mode using the IRS low-resolution modules, 
Short-Low (SL) and Long-Low (LL), which span wave- 
lengths from 5-14 /zm and 14-38 /im, respectively, with 
a resolution X/5X ~90. 

Details on the observational tec hniques and genera l 
dat a, reduction steps can be found in lFurlan et al.l (2006) 
and lWatson et al.r(|2009[ ). We provide a brief summary. 
Each object was observed twice along the slit, at a third 
of the slit length from the top and bottom edges of the 
slit. Basic calibrated data (BCD) with pipeline version 
S18.7 were obtained from the Spitzer Science Center. 
With the BCDs, we extracted and calibrated the spectra 
using the SMART package (|Higdon et al.l 120041 ) . Bad 
and rogue pixels were corrected by interpolating from 
neighboring pixels. 

Most o f the data were sky subtr acted using optimal ex- 
traction dLcboutcil ler et aLll2010[) . The exceptions were 
LRLL 37, LRLL 55, and LRLL 68. In LRLL 37, there 
was an artificial structure in the 5-8 /jm region which 
was removed by performing off-nod sky subtraction. A 
similar structure was seen in LRLL 55 and LRLL 68, but 
due to the high background in the area, optimal extrac- 
tion was necessary for the LL order. Therefore, the final 
spectra for LRLL 55 and LRLL 68 are a combination of 
the off-nod sky-subtracted SL spectra and the optimally 
extracted LL spectra. 

To flux calibrate the observations we used a spectrum 
of a Lac (Al V). We performed a nod-by-nod division of 

11 http:/ /www. astronomy.ohio-state.edu/~jdeast/4k 
/proc4k.pro 

12 http://web.mit.edu/~burles/www/MIKE/ 



the target spectra and the a Lac spectrum a nd then mul- 
tiplie d the result by a template spectrum ([Cohen et al.l 
I2003fh The final spectrum was produced by averaging the 
calibrated spectra from the two nods. Our spectropho- 
tometric accuracy is 2-5% estimated from half the dif- 
ference between the nodded observations, which is con- 
firmed by comparison with IRAC and MIPS photometry. 
We note that there are artifacts in the spectra of LRLL 68 
and LRLL 133 beyond ^30 /im. For clarity, we manually 
trim the spectra to exclude these regions. The final spec- 
tra used in this study are shown in Figures El 02 and [4] 

2.4. Millimeter Flux Densities 

We observed LRLL 21, LRLL 31, LRLL 67, LRLL 68, 
and LRLL 72 in IC 348 with the Sub-millimeter Array 
(SMA) on November 12, 2008. We used the Compact 
Configuration with six of the 6 meter diameter anten- 
nas at 345 GHz (860 /im) with a full correlator band- 
width of 2 GHz. Calibration of the visibility phases 
and amplitudes was achieved with observations of the 
quasars 3C 111 and 3C 84, typically at intervals of 20 
minutes. Observations of Uranus provided the absolute 
scale for the flux density calibration. The data were 
calibrated using the MIR software package We de- 
tected LRLL 31 and LRLL 67 with flux densities of 
0.062±0.006 Jy and 0.025±0.011 Jy, respectively. We did 
not detect LRLL 21, LRLL 68, or LRLL 72 and measure 
a 3<7 upper limit of 0.015 Jy for these objects. 

3. ANALYSIS 

Here we model the SEDs of the targets in our sample. 
First we collect the stellar properties of our objects, ei- 
ther by adopting literature values or deriving our own in 
Section [3. II These stellar properties are important input 
parameters for our physically motivated models which 
we discuss in Section [3T21 In Section |3~31 we discuss the 
results of our SED model fitting, as well as the degenera- 
cies that exist given that we lack millimeter observations 
for many of our targets. 

3.1. Stellar Properties 

Stellar parameters are listed in Ta ble [3] M« was de- 
rived from the HR diagram and the iSiess et al.l (|2000D 
evolutionary tracks using T* and L*. Stell ar tem- 
peratures are from iKenyon fc Hartmannl (|1995| ) , based 
upon the spectral types adopted for the targets in 
NGC 2068 and IC 348 (iFlahertv fc Muzerollel [20081: 
iLuhman et al.l 120031 respectively). Luminosities are 
cal culated with dereddened J-b and photometry follow- 
ing [Kenyoir&|Tla7tma^^ of 
400 pc for NGC 2068 dFlahertv fc Muzerollel l2008h and 
315 pc for IC 348 ([Luhman et al.ll2003l ). R* is calculated 
using the derived luminosity and adopted temperature. 
The derivation of the mass accretion rates and accretion 
luminosities are discussed in Section [3. 1.11 

All photometry for our NGC 2068 objects is taken 
from lFlahertv fc Muzerollel ([2008I L This includes BVRI, 
2MASS JHK, IRAC and M IPS data. We note that 
iFlahertv fc Muzero lle (2008) present photometry from 
the SDSS surve y, which we c onver t to Johnson- Cousins 
BVRI following IJordi et al.l (l2006lh IRAC and MIPS 

13 http:/ /www. cfa.harvard.edu/~cqi/mircook. html 
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photometry for IC 348 comes from lLada et all (12006 
and 2MASS JHK data is from Skrutskie et all (2006 



UBVRI photometry for our IC 348 targets comes mainly 
from this work, but is supplemented by values in the lit- 
erature. All UBVRI photometry for LRLL 31, LRLL 55, 
LRLL 67, LRLL 68, LRLL 72, and LRLL 133 are solely 
from this work. UVR data for LRLL 21 and LRLL 37 
are from this work, b ut we use I-band magnitudes from 
iLuhman et all (2003) for b oth a nd a B-band magnitude 
for L RLL 21 from |Hgrbij {1998). We use R and I data 
from lHerbid (fl998h f or LR LL 2. BVRI data for LRLL 6 
is als o fromlHerbis) (1 19981) . All L-band magnitudes are 
from lHaisch et all (|200Tf K 

Extinctions were measured by comparing V-R, V- 
I, R-I, and I-J colors to photospheri c colors from 
IKenvon fc Hartmannl (|1995D . We used the lMathisI (|1990f ) 
extinct ion law for obje cts with Ay< 3. For Ay> 3, we 
use the lMcClurd (|2009D extinction law. We adopt Ry=5 
which is more approp riate for the denser regions stud- 
ied here (|Mathislll990l) . In most cases, extinctions based 
on I-J colors gave the best fit. Since we have no I-band 
data for FM 581 we adopt an extinction based on V-R; 
LRLL 2 and LRLL 6 have no VRI photometry in the 
literature and so we adopt the extinction measured by 
ILuhman et all (|2003l ). All extinctions used in this work 
are listed in Table |3] We note that most of our extinc- 
tions are similar to those in the lite rature (Luhma n et all 
120031: iFlahertv fc Muzerolld 120081) . In cases where we 
found differences, we chose to rely on our measurements 
since they are based on I-J colors which have recently 
been shown to be the l east affected by exc ess emission 
at shorter wavelengths (|Fischer et al.l |20TTI McClure et 
al., in preparation). For early- type stars, the peak of the 
stellar emission would be at shorter wavelengths. How- 
ever, the early-type star in our sample (LRLL 2) has bad 
photometry so we do not explore this point further. 

3.1.1. Accretion Rates 

During the classical TTS phase of stellar evolution, 
young objects accrete material from the disk onto the 
star via magnetospheric accretion (jUchida fc Shibatal 
1984). The infalling gas impacts the stellar surface at ap- 
proximately the free fall velo city creating a shock whic h 
heats the gas to ~ 1 MK (|Calvet fc Gullbringj Il998l ). 
The shock emission is reprocessed in the accretion col- 
umn and the observed spe ctrum peaks in the ultraviolet 
(jCalvet fc Gullbrmdll998l) . The best estimate of the ac- 
cretion rate is found by measuring the total luminosity 
emitted in the accretion shock, i.e. the accretion lumi- 
nosity. While ultraviolet emission is difficult to observe 
from ground-based observatories, there are many tracers 
of the accretion lum inosity at longer wavelengths (see 
IRigliaco et all (|201 If ) for a comprehensive list). A few of 
those tracers include excess emission observed in U-band 
photometry, emission in the NIR Ca II triplet lines, and 
the Ha line profile. 

Second to the ultraviolet excess, U-band excesses are 
the best measure of the flux produced in the shock and 
have been shown to correlate w ith the total shock ex- 
cess (jCalvet fc Gullbringj 119981 ) . However, at low M 
emission at U may be dominated by chromospheri c 
emission (jHoudebine et all 119961: iFranchini et~all fl998): 
this chromospheric excess can confuse d eterminations of 
the accretion rate (|Inglebv et al.l 120111) . While possi- 



ble from the ground, U-band observations are still dif- 
ficult to obtain, especially when extinction towards the 
source is relatively high as in the case of our sample, 
all with Ay > 1. Observations of the Ca II near- 
infrared triplet are easily obtained from the ground, 
even for high Ay sources, and the flux in the 8542 
A line also correlate s with the accretion luminosity 
(|Muzerolle et al.l [19981 Ca II is observed m emission m 
accreting sources but is unreliable at low accretion rates, 
when the chrom ospheric emission rivals that from accre- 
tion in strength (lYang et"aT] [2007: Batalha & Basij HT993l : 
llnglebv fc et alll201lD . Ha is commonly used as a tracer 
of accretion, both by measuring the equiv alent width and 
the ve locity width of the line in the wings (IWhite fc Basri 
2003t IBarrado v Navascues fc Ma rtin 20031: INatta et al 



20041 ). Models of magnetospheric accretion have repro- 



duced the observed velocities in Ha, tracing material 
traveli ng at several hundred km s -1 near the accretion 
shock (|Lima et all 1201 CI: IMuzerolle et al.l 1 2001 . 

The mass accretion rates adopted for our sample are 
listed in Table [3] We use d U-band photometry from this 
work and the relation in lGullbring et all ([ 19981 ) to mea- 
sure mass accretion rates for LRLL 37 and LRLL 68 and 
an upper limit for LRLL 1330 For FM 177, FM 281, 
FM 515, FM 581, FM 618, LRLL 21, LRLL 67, and 
LRLL 72, we used high resolution echelle spectra ob- 
tained with MIKE to measure the width of the Ha line in 
the wings (at 10% of the maximum flux; Figure [lj. We 
then compared this to the relation between line width 
and M in INatta et all (|2004D to obtain mass accretion 
rate estimates for these eight sources. While our MIKE 
spectra covered both Ha and the Ca II triplet, Ha pro- 
vided a more accurate estimate of M due to confusion 
with chromospheric Ca II emission at the levels of accre- 
tion found in these sources. Given the chromospheric ap- 
pearance of the Ha profile of LRLL 72, its mass accretion 
rate should be taken as an upper limit. For LRLL 31 we 
adopted a mass accretion rate from the literature. We do 
not have mass accretion rate measurements for LRLL 2 
or LRLL 6. 

For NGC 206 8, we compared our derived accretion 
rates to those in Flahe rty fc Muzerolld (|2008l) who cal- 
culated the amount of excess continuum emission nec- 
essary to produce the observed veiling of the photo- 
spheric lines. Within a factor of 2-3, the normal er- 
ror in M estimations, both calculations of the accretion 
rate agree, with a few exceptions. When comparing the 
Ha line widths at 10% we find that our MIKE line pro- 
file of FM 281 is ~ 180 kms" 1 narrow er than when 
observed by IFlahertv fc MuzeroTIe] (|2008f ). In addition, 
our observation of FM 177 is con sistent with an accret- 
ing so urce, while when observed bv lFlahertv fc Muzerolld 
(2008) its Ha profile was consistent with that of chromo- 
spheric emission. Variability is known to occur in T Tauri 
stars so the decrease in Ha and M are not unexpected 
(|Codv fc H illenbrand 2010). For these objects we chose 
to adopt the accretion rate obtained from our MIKE ob- 
servations. The biggest uncertainty in calculating accre- 



We measured an upper limit for LRLL 55 of 4xl0" 6 M Q 
yr — 1 using U-band photometry. However, this very high upper 
limit does not provide useful constraints for the purposes of this 
paper and so we do not comment on it further. 
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tion rates using veiling is the choice of bolometric cor- 
rection, which can vary in value by a factor of 10 de- 
pend ing on which analysis is used (|White &: H illenbrand 
120041 ) and the spectrum of the excess emissi on which veils 
the p hotospheric lines can be complicated ((Fischer et al.1 

1ml). 

3.2. Disk Model 

We try to reproduce the SEDs presented in Fig- 
ures [2l [3l & H] usin g the i r radiated accretion disk m odels 
of iD'Alessio et"afl (fl99l [19991 [200l [2001 I2001I . We 
point t he reader to those pap ers for details of the model 
and to lEspaillat et"ail (|2010l) for a summary of how we 
apply the model to the SEDs of transitional and pre- 
transitional disks. Here we provide a brief review of the 
salient points of the above works. 

When we refer to a "full disk model" we mean a disk 
model composed of an irradiated accretion disk and a 
frontally illuminated wall at the inner edge of the disk 
which is located at the dust sublimation radius. The 
inner wall dominates the emission in the NIR, the wall 
and disk both contribute to the MIR emission, and the 
disk dominates the emission at longer wavelengths. Com- 
pared to a full disk model, a pre-transitional disk model 
has a gap within the disk. In this case, we include a 
frontally illuminated wall at the dust sublimation ra- 
dius and another wall at the gap's outer edge. For this 
outer wall we include the shadow cast by the inner wall 
(jEspaillat et al.ll201Q[ ). We do not include an inner irra- 
diated accretion disk behind the inner wall since previous 
work has shown that the inner wa ll dominates the emis - 
sion at these shorter wavelengths (jEspaillat et al.ll2010[) . 
Behind the outer wall we include an irradiated accretion 
disk in cases where we have millimeter data, which is 
necessary to constrain the outer disk. The inner wall 
dominates the NIR emission while the outer wall dom- 
inates the emission from ~20-30 fim. The outer disk 
dominates the emission beyond ~40 /im. A transitional 
disk model is very similar to that of a pre-transitional 
disk model except that we do not include an inner wall 
at the dust sublimation radius. In some instances, we in- 
clude a small amount of optically thin dust in the inner 
hole or gap in transitional and pre-transitional disks to 
reproduce the 10 /jm silicate emission feature. We cal- 
culate th e emission from thi s optically thin dust region 
following ICalvet et al.l (|200l . 

3.2.1. Disk Properties 

Table [?] lists the model-derived properties of our sam- 
ple. The heights of the inner and outer walls (z wa u) and 
the maximum grain sizes (a maa; ) are adjusted to fit the 
SED. T wa u is the temperature at the surface of the opti- 
cally thin wall atmosphere. The temperature of the in- 
ner wall of full disks and pre-transitional disks (T', )n H ) is 
held fixed at 1400 K (except for FM 515, see Section GO]) 
which is the typical temper ature o f dust at the sublima- 
tion radius (jMuzerolle et al.ll2003f) . The temperature of 
the outer wall (T^aii) m transitional and pre-transitional 
disks is varied to fit the SED, particularly the IRS spec- 
trum. The radius of the wall (R wa ii ) is d erived using the 
best fi tting T wa u following Equation 2 in Espaillat et al.l 
(|2010h . 

Previously, we have seen that in transitional and pre- 
transitional disks, the IRS spectrum is dominated by the 



outer wall while the outer disk dominates the millimeter 
emission (Ell). Since the majority of the emission seen 
by IRS is from the outer wall, the IRS spectrum is a good 
constraint of the hole/gap size and in most cases SED- 
derived hole/gap sizes are in reasonable agreement with 
those obtained with millimeter imaging ([Andrews et al.l 
1201 If ). For objects in this work with no millimeter data, 
we do not include an outer disk given that its contribu- 
tion to the IRS SED is expected to be small. We did 
include an outer disk for the two objects in the sample 
for which we have millimeter data: the transitional disk 
LRLL 67 and the pre-transitional disk LRLL 31. We also 
included an outer disk when modeling each of our full 
disk targets since additional emission from an outer disk 
is necessary to reproduce the observed MIR emission. 
The parameters of the outer disk which are varied to fit 
the SED are the viscosity parameter (a) and the settling 
parameter (e; see Section [3. 2. 2|) . One can interpret vary- 
ing a as fitting for the disk mass since Mdisk oc M/a (see 
Equation 38 in ID'Alessio et al.l ll998). We note that the 
mass accretion rate onto the star does not necessarily re- 
flect the mass transport across the outer disk, especially 
in the case of TD and PTD where the mass accretion 
rate onto the star is likely an underestimate of the mass 
transport across the outer disk (see Section 14. 4[) . We 
also do not expect that the mass accretion rate is con- 
stant throughout the disk. However, for simplicity, here 
we assume that the mass accretion rate measured onto 
the star is representative of the disk's accretion rate. 

We will discuss how the lack of millimeter constraints 
leads t o de generacies in our outer disk model fitting in 
Section [3~3l There we also discuss the effect the adopted 
disk inclination and outer radius have on the simulated 
SED. We assume that the inclination of the disk is 60° 
for all of our objects and that they have an outer disk 
radius of 300 AU (except FM 581, see Section [373]) . 

3.2.2. Dust Properties 

The opacity of the disk, and hence the temperature 
structure and resulting emission, is controlled by dust. 
The dust opacity depends on the composition of the dust 
assumed. It also depends on changes in the dust due to 
grain growth and settling. Grains grow through colli- 
sional coagulation and settle to the disk midplane due to 
gravity. Since in this work we include models of several 
full disks, here we review the effect that the dust prop- 
erties have on t h e disk structure in more detail following 
ID'Alessio etldl (|2006l ). 

Dust Settling. The settling of dust has important, 
and often overlooked, effects on the disk's density- 
temperature distribution and emission. When there is 
some degree of settling, the dust-to-gas mass ratio of 
grains in the disk atmosphere decreases with respect to 
the standard value (i.e. the diffuse interstellar medium) . 
This has several effects: (1) it decreases the opacity of 
the upper layers; this allows the impinging external ra- 
diation to penetrate deeper into the disk, decreasing the 
height of the irradiation surfaced and making it geo - 
metrically flatter (see Figure 3 in ID'Alessio et al.ll2006f h 
which in turn decreases the fraction of the irradiation 

15 The height of the disk irradiation surface, z s , is defined by the 
region where t s , the radial optical depth to the stellar radiation, is 
~ 1. 
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flux intercepted by this surface^ decreasing the contin- 
uum flux emerging from the disk, (2) since most of the 
external radiation is deposited at the irradiation surface, 
lowering it changes the temperature-density structure of 
the atmospheric layers, where the temperature inversion 
occurs (also called the super- heated layers), modifying 
their contribution to the SED, and finally, (3) it changes 
the emissivity of the disk interior; in this region the dust- 
to-gas mass ratio of the grains increases given that the 
grains removed by depletion from the upper layers are 
now located deeper in the disk. 



Some of the above-mentioned effects can be ac- 
counted for by arbitrarily changi ng the disk surface 



height as a function of radius (e.g. Mivakc & Nakagawa 


1995 


:lCurrie & Sicilia-AguilaHl201lMSicilia-Aguilar et al. 


2011 


), and this will probably give a reasonable estimate 



of the continuum SED of the disk. However, the con- 
tribution of the upper layers to the SED or the role of 
the deeper layers in millimeter images and emergent flux, 
would not be consistent for this simple approach to set- 
tling. On the other hand, taking into account the de- 
tailed physics of settling (e.g. IWe idcnschill ing et al .111 9971 : 
IDullemond &: D ominik 2004) is complex and simulations 
show that disks should be completely settled within ^10 6 
years, in contradiction with observations, reflecting that 
we are missing some processes that can keep some small 
grains in t he upper layers for longer t imescales (e.g. 
turbu lence; IDullemond fc D ominik 2005; Birnstic l et al.l 
12011 . 

In our SED m odeling we have adopt ed a different ap- 
proach following ID'Alessio et al.l ([2006D by parameteriz- 
ing settling as a depletion of dust in the upper layers, 
with a corresponding increment of the dust-to-gas mass 
ratio near the midplane. The maximum grain sizes in the 
disk atmosphere and interior are allowed to change, re- 
flecting the possibility of grain growth. We can also vary 
the height in the disk that separates the atmosphere from 
the interior as well as the degree of settling. The amount 
of settling is parameterized by e = Catm/Cstd, (i-e., the 
dust-to-gas mass ratio of the disk atmosphere divided by 
the standard value). The main point of this approach is 
that the same grains that determine the height and shape 
of the irradiation surface and the amount of intercepted 
external flux, are the ones that are emitting in the mid- 
IR silicate bands, and their emissivity and temperature 
distribution are consistent with their properties. Also, 
the grains near the midplane which are responsible for 
the mm emergent intensity have a dust-to-gas mass ra- 
tio related to the properties of the atmospheric grains. 
The advantage of such an approach is that, in principle, 
observations can be used to constrain the grains' com- 
position, size and spatial distribution, and this can be 
related to models of the detailed dust evolution in disks. 
However, to really fulfill this goal, we need observations 
that cover a wide range of wavelengths with high reso- 
lution. Given our present observations, we have chosen 
to adopt a radially constant e and to assume that the 

16 As stellar radiation enters the disk, it does so at an angle 
to the normal of the disk surface (6>o = cos ~ 1 fio). A fraction of 
the stellar radiation is scattered and the st ellar radiatio n captu red 
by the disk is ~ (<rT*/ir){R* /R) 2 fi (see ICalvet et al.l (fl99lT) for 
further discussion). Therefore, if the disk is more flared, fiQ is 
larger and more stellar irradiation will be intercepted by the disk 
and it will be hotter and emit more radiation. 



interior grains are concentrated very close to the mid- 
plane (at z < 0.1 H). These assumption s will not affect 
the mid-IR SED (|D'Alessio et al.l 120061 Qi et al. 2011) 
and we avoid introducing new sets of free parameters to 
the problem, retaining the important physical properties 
of settling. 

Dust Grain Growth. In this work we also change the 
maximum grain size in the disk. The models assume 
spherical grains with a distribution of a~ p where a is 
the grain radius be tween a m „ and a, max and p is 3.5 
(jMathis et alJfT9~77h . A mixture with a smaller a max has 
a larger opacity at shorter wavelengths than a mixture 
with a larger & max . Since the height of the disk surface, 
z s , is defined by the region where t s ~ 1 small grains will 
reach this limit higher in the disk relative to big grains. 
Therefore, disks with a small B, max are more flared than 
disks with a large a max for the same dust-to-gas mass 
ratio. One difference between increasing the settling and 
increasing the grain size is that with settling, small grains 
remain in the upper disk layers and so we still see silicate 
emission while with grain growth in the disk atmosphere, 
the silicate emission disappears since larger grains do not 
have this feature in their opacity. In the walls and the 
outer disk, a m i n is held fixed at 0.005 /im while a max is 
varied between 0.25 /xm and 10 /zm to achieve the best 
fit to the silicate emission features. In the outer disk, 
there are two dust grain size distributions as mentioned 
above. In the disk interior the maximum grain size is 
1 mm (jD'Alessio et al.l 120061 ) . The maximum grain size 
of the disk atmosphere is adjusted as noted earlier. 

Dust Composition. The composition of dust used in 
the disk model imp acts the resulting SEP and derived 
disk properties (see lEspaillat et al.l (|2010D for a discus- 
sion). We follow Ell and perform a detailed dust compo- 
sition fit for the silicates seen in the IRS spectra including 
olivines, pyroxenes, forsterite, enstatite, and silica. We 
list the derived silicate mass fractions in Tables [5] and [7] of 
the Appendix. In addition to silicates, we also included 
organ ics, troilite, and water ice following Espaillat et al.l 
(2010) and Ell. We note that only silicates exist at 
the high temperatures at which the inner wall is located. 
In transitional and pre-transitional objects where we in- 
clude optically thin dust within the hole, the silicate dust 
composition and abundances are listed in Tables [5] and [3] 
of the Appendix, respectively. 

3.3. SED Modeling 

In this work we present the first detailed modeling of 
disks with IRS spectra in NGC 2068 and IC 348. We find 
that all the objects are reasonably reproduced with tran- 
sitional, pre-transitional, or full disk models. It is not the 
goal of this paper to find a unique fit to the SED. To ar- 
rive at a unique fit, one would ideally have a finely sam- 
pled, multi- wavelength SED as well as spatially resolved 
data at multiple wavelengths. Finely sampled data on 
the time domain would also be necessary since the emis- 
sion o f TTS is known to be variable (e.g. Espaillat et al. 
2011). Given that this situation is not currently achiev- 
able, here we focus on finding a fit that is consistent 
with the observations presented in this work. Our as- 
sumptions of the disk structure are an oversimplification. 
Recent hydrodynamical simulations show that the inner 
regions of tr ansitional and pre-transitional disks should 
be complex (|Zhu et al.ll20lH iDodson-Robinson fc Salvkl 
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1201 lfl . However, in the absence of data capable of con- 
firming these simulations, we proceed with our simple 
model. We discuss additional assumptions and how they 
play into the degeneracies of our modeling in Section f4.2l 
We present details of the derived dust composition in the 
Appendix. 

3.3.1. Results 

We find a large range of hole and gap sizes for our 
transitional and pre-transitional disks. Our TD have 
holes spanning 4 to 49 AU (Figure EJ Table |4}. Our 
pre-transitional disk targets have gaps ra ngin g from 5- 
45 AU (Table |1 and Figures B H and [7]). The three 
objects easily identified by dips in the SED (FM 515, 
FM 618, LRLL 31) have gap sizes of 11-45 AU. We note 
that here we classify LRLL 21 as a PTD even though 
its NIR emission is weaker than the other PTD in our 
sam ple and it resembles t he emission expected from a 
TD. iFlahertv etail ((20T2T ) find that LRLL 21 has sig- 
nificant NIR emission in more recent IRS observations, 
pointing to strong intrinsic variab ility in the inner dis k 
linked to changes in the inner wall (jEspaillat et al.ll2011[ ). 
Therefore, here we classify LRLL 21 as a PTD. Another 
PTD in our sample that is not obvious based on its SED 
alone is LRLL 37 which has the smallest gap size in our 
sample (5 AU; Figure [S]). It is not possible to fit the 
IRS data of LRLL 37 with a full disk model, even within 
the uncertainties of the observations. In particular, we 
could not fit the strong 10 /im silicate emission with our 
full disk model. This could be a sign that LRLL 37 is 
a pre-transitional disk with a small gap that contains 
some small optically thin dust, reminiscent of RY Tau 
(Ell) and we will return to this point in Section |4~T1 We 
only have millimeter fluxes for two objects in our sample, 
LRLL 31 and LRLL 67. For these disks we derive a disk 
mass of 0.06 M Q for each (Figure [7]). The best fitting e 
and a for LRLL 31 arc 0.001 and 0.005, respectively. For 
LRLL 67, e=0.001 and a=4xl0~ 5 . 

Most of the transitional and pre-transitional disks have 
small optically thin dust within the inner 1 AU of the 
hole or gap. The exceptions are the transitional disk 
LRLL 72, where we find the 10 /im silicate emission can 
be produced by the optically thin atmosphere of the outer 
wall, and the pre-transitional disk LRLL 31, where the 
10 //m silicate emission comes from the inner wall's at- 
mosphere and optically thin dust within the gap is not 
necessary to fit the observations. The mass of dust and 
sizes of the grains in this region are given in the Ap- 
pendix. 

For FM 581, LRLL 2, LRLL 6, LRLL 55, and LRLL 68 
we can fit the SED reasonably well using full disk mod- 
els (Figure HI). Unlike LRLL 31 and LRLL 67 above, we 
do not have millimeter detections for FM 581, LRLL 2, 
LRLL 6, LRLL 55, and LRLL 68 and so we cannot con- 
strain the mass of the disk. Therefore, the models pre- 
sented here are more uncertain, but we show them to 
illustrate that a disk model with dust settling and no 
holes or gaps in the disk can reproduce the observed 

17 As mentioned in S ection 13.21 we do not i nclude an inner disk 
behind the inner wall. lEspaillat et ail 120101) find th a t the inner 
wall dominates the NIR emission of PTD. IPott et~aTI iwWi ) con- 
firm that the inner disk in PTD is small. High resolution imaging 
is needed to further constrain the gap sizes of the objects presented 
in this work. 



SEDs. We do have an upper limit for the millimeter 
flux of LRLL 68 and we use this object as an example 
to discuss the degeneracies inherent in the modeling pre- 
sented, mainly due to lack of millimeter detections, in 
Section 13.3.21 To briefly summarize, disk models with 
the same e-to-a ratio will produce very similar emission 
in the IR but substantially different emission in the mil- 
limeter. Therefore, millimeter data is crucial to disen- 
tangle this degeneracy and the disk parameters in this 
work should only be taken as indicative of a model that 
can reproduce the observed SEDs. 

With the above degeneracies in mind, we limited our 
parameter search and set e=0.001, changing only a until 
we achieved a good fit to the SED. We could have also 
set a to a certain value and fit for e instead, however, as 
mentioned above and as discussed in Section 13.3.21 e/a 
is the most relevant result. For LRLL 2, LRLL 55, and 
LRLL 68 we find a=0.06, 0.004, 0.006, respectively. For 
LRLL 6 using an e of 0.001 required a>0.1 which would 
lead to a viscous timescale shor ter than the lifetime of 
the disk (|Hartmann et al.l IT998;) . so instead we set e= 
0.0001; the best-fitting a in this case was 0.1. To fit 
the very steep downward slope of FM 581, we needed 
to significantly truncate the outer disk radius, down to 
0.6 AUE3 We fit FM 581 with e=0.001 and a=0.00006. 

3.3.2. Model Degeneracies and Millimeter Constraints 

Here we explore the degeneracies introduced into the 
modeling presented in this paper due to the lack of mil- 
limeter data. Rather than do this for each object, we 
selected LRLL 68 for this test since it has an upper limit 
to its millimeter flux at 860 /im from the SMA observa- 
tions reported in this paper. 

First, disk models (around the same star) with the 
same e-to-a ratio will have SEDs with similar emission in 
the IR. This is because disks with equal e/a have similar 
disk surfaces. The disk surface is defined as the point 
in the upper disk layers where the radial optical depth 
(which depends on the product of the opacity and col- 
umn density) to the stellar radiation reaches one. e de- 
termines the abundance of small dust (i.e. the opacity of 
the upper disk layers) while a affects the surface density. 
Therefore, in disks with equal e/a the same fraction of 
stellar flux will be intercepted by the disk. Since the IR 
is dominated by the upper layers of the inner disk, their 
emergent intensity in the IR will be similar. 

While the IR emission is similar, the emission seen 
in the millimeter will be different. A disk with an 
a— viscosity has a mass surface density given by S « 
Mflk/a < c s >, where f2& is the Keplerian angular ve- 
locity and < c s > is the sound speed. This implies that 
the disk mass, for a given disk radius and similar outer 
disk temperature distribution, would be proportional to 
a . Therefore, a disk with a smaller a will have a larger 
disk mass. In the millimeter, we are more sensitive to the 

18 T he SEP of FM 581 r esembles that of the 5-10 AU binary 
SR 20 (IMcClure et al.ll200B . SR 20's disk is outwardly truncated 
at 0.4 AU, too far to be attributable to the known companion, and 
so this trunca tion would have to b e due to an unseen companion 
at ~l-2 AU HMcClure et al.ll2008T) . Likewise, it seems that the 
most viable mechanism to truncate the disk of FM 581 to such 
small radii would also be a companion. However, FM 581 still has 
a substantial accretion rate. Millimeter observations are necessary 
to constrain the mass and size of this disk. 
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big grains in the midplane of the disk, where most of the 
disk's mass is stored, and so disks with small a and a 
higher column density will have more millimeter emis- 
sion. This means that models with similar ecos(i)M/a 
would have similar IR SEDs, but different millimeter 
SEDs. Millimeter data is necessary to disentangle this 
degeneracy. 

Because of the above, we can find a best-fit e/a to the 
IR emission if we hold i and M constant. However, each 
model will have different emission in the millimeter and 
since we have no millimeter data, we cannot claim that 
a particular combination of e and a is better than an- 
other. For the models shown in Figure [9] (Models 1, 2, 
& 3 in Table [5]) we hold the mass accretion rate, stel- 
lar parameters, inclination (i), and outer disk radius Rd 
fixed and change only e and a. We set e to 0.0001 , 0.001 , 
0.01, and 0.1, the values given in iD'Alessio etldl ([2006), 
and fit the SED by changing a. We do not discuss cases 
where a>0 .1 since these disks wo uld have short viscous 
timescales ( Hartmann et al.lll998l ). For LRLL 68, we find 
that the best- fitting e/a is 0.2. 

We also look at how varying the outer disk radius and 
inclination affect the simulated SED. Changing the ra- 
dius changes the mass of the disk but the millimeter 
emission does not change significantly (Figure fTO)) . This 
is because the mass depends on the disk radius (M^ = 

f R d T 1 2ttR(1R). However, the column density of the an- 
nuli contributing to the mm flux remains the same. This 
highlights that disk sizes cannot be firmly constrained 
without resolved imaging. We find that changing the 
inclination angle while holding other parameters fixed 
changes the IR emission, but does not significantly alter 
the millimeter emission (Figure [TTj) . The disk is mostly 
optically thin at millimeter wavelengths so we can see 
through the disk at any inclination. Therefore, the mil- 
limeter emission does not change with inclination. On 
the other hand, the IR emission does change; as the in- 
clination decreases, we see more of the inner disk, which 
mainly dominates the IR emission. We note that the 
near-IR emission also depends on the shape of the wall. 
We assume the wall is vertical, therefore it will not con- 
tribute to the SED at 0° or 90° and will p roduce the 
most emission at 60° (jDullemond et aLll2001l ). If the wall 
is curved, we would expect to see more emission at lower 
inclinations (|Isella fc Nattall2005l) . 

4. DISCUSSION 

4.1. Limitations of Observations 

As discussed bv lEspaillat et al.l ([2010D . the sizes of gaps 
and holes we can detect in the disk are limited by broad- 
band SEDs. Given that the majority of the emission at 
~10 /im in a typic al disk traces the dust within the inne r 
1 AU of the disk (jD'Alessio et al.l[2006l lEspaillat][2009l ). 
the Spitzer IRS instrument will be most sensitive to clear- 
ings in which much of the dust located at radii <1 AU 
has been removed. Because of this, IRS is more effective 
in picking out disk holes where dust at small radii has 
been removed. However, IRS cannot easily detect gaps 
whose inner boundary is outside of ~1 AU. For exam- 
ple, a disk with a gap ranging from 5 - 10 AU will be 
diffic ult to distinguish from a full disk (Espaillat et al. 
2010). The gaps currently inferred solely from SEDs, 
which have been modeled and imaged in the millimeter, 



are typically quite large. This reflects an observational 
bias towards picking out pre-transitional disks with large 
gaps since their mid-infrared deficits will be more obvious 
in Spitzer spectra. Smaller gaps will not have as obvious 
of a deficit and will be difficult to detect. Broad-band 
colors from IRAC and MIPS have their limitations as 
well. They are useful for picking out TD, but it is diffi- 
cult to distinguish the NIR colors of a PTD from a full 
disk. 

LRLL 37 may be a case of a disk with a small 5 AU 
gap. The dip in the SED is not obvious, but the strong 
silicate emission seen in this object is reminiscent of RY 
Tau. RY Tau has a large cavity in its disk based on 
millimeter imaging ([Isella et al.l (201(1 and its SED was 
modeled with a gap of ~20 AU (Ell). There are many 
other disks that ex hibit strong 10 /im silicate emission 
(jFurlan et al.l [20 09) and perhaps this is a hint pointing 
to small gaps in disks. That said, we cannot exclude 
the presence of small gaps in what we have labeled full 
disks in this paper, or for that matter any full disks in 
general. High-resolution imaging is crucial to investigate 
this further and current TD fractions should be taken as 
lower limits. 

4.2. Limitations of Models 

The underlying physical structure inferred from the 
SED is dependent on the model one uses. For ex- 
ample, there are five disks in Taurus with decreas- 
ing emission at IR wavelengths, much like the disks 
studied in this work, that hav e seemingly contrast - 
ing interpretations presented by ILuhman et all (2010) 
and iCurrie fc Sicilia-Aguilarl ([201 ll ). The main differ- 
ence between the two papers are the models adopted. 
ICurrie fc Sicil ia-Aauil an d20ll use the model grid pre- 
sented in iRobitaille et all (j'2007) . The authors account 
for the observed SEDs by decreasing the mass of models 
with well-mixed gas and dust in the disk (i.e. a disk with 
no settling) to the point where it be comes entirely opti- 
cally thin. The ILuhman et ail (|2010f) re sults are based o n 
synthetic colors from the model grid of lEspaillatl ((2009) . 
Those models are the s ame in this wo r k, whe re dust set- 
tling is incorporated. ILuhman et al.1 (|2010f ) can repro- 
duce the observed colors with disks that have dust set- 
tling. To illustrate that a settled disk can reproduce the 
broad-band SED as well, in the Appendix we model ZZ 
Tau, one of the five disks in question. Therefore, it be- 
comes clear that the different interpretation between the 
two groups is biased by the models adopted. The most 
one can say is that both a well-mixed low-mass disk and 
a settled disk can reproduce the observations. In the for- 
mer case the disk is optically thin to its own radiation 
at all radii; for the latter case the innermost disk is still 
o ptically thick to its own radi ation. 

iSicilia-Aguilar et al.1 ([201 ll ) also independently find 
that they need to incorporate settling to reproduce the 
SED of objects with decreasing IR emission. We note 
that the implementation of settling in that wo r k and 
our work are very different. ISicilia-Aguilar et al.l ([20111 ) 
simulate settling in their modeling by lowering the disk 
surface, but its shape remains the same. In addition, 
the dust-to-gas mass ratio is held fixed throughout the 
disk (i.e. there is no dust depletion in the upper disk 
layers). Since the surface is still flared and the opacity 
remains the same, the disk is hotter than it would be if 
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the disk was geometrically flatter and the opacity was 
lower. Therefore, such a disk will produce more emission 
and there will be an inherent bias towards decreasing the 
disk mass in order to reproduce lower observed fluxes. In 
the models used in this work we deplete the small grains 
in the upper atmosphere of the disk and self-consistently 
calculate the disk height and shape. This is an iterative 
process given that the height of the disk irradiation sur- 
face dictates how much energy is captured by the disk 
and this depends on the opacity set by the dust proper- 
ties and the density of the disk, which in turn depends on 
the temperature through the scale height. We note that 
it is still possible to have a disk which is both settled 
and low-mass. Our point is that in order to constrain 
the dis k mas s and avoid the degeneracies discussed in 
Section 13.3.21 millimeter data is necessary and the full 
effects of dust settling need to be taken into account. 

This leads to another issue that is quite model de- 
pendent: the disk mass. Not only does the disk mass 
depend on the opacity one assumes, it also depends on 
the surface density and temperature radial profiles. For 
example, our mass d eterminations (i.e. Ell) ar e con- 
sistently higher than lAndrews fc Williams (2005) since 
our opacities are ~3 times lower, we assume the outer 
disk radius is larger, and we use a self-consistent surface 
density and temperature instead of power-law approxi- 
mations. Therefore, masses obtained by different models 
cannot be meaningfully compared. 

4.3. Disk Structure & Semantics 

There are many terms in the literature aiming to cat- 
egorize SEDs of TTS surrounded by disks. At times this 
leads to confusion. For example, what some researchers 
call a transitional disk others would call an anemic disk or 
an evolved disk or a homologously depleted disk. The ef- 
fect of not consistently applying the term "transitional" 
in the literature is seen when looking at TD fractions 
in the literature. Some call "transitional" any objects 
whose SED does not resemble the median SED of Tau- 
rus. Others use the more restrictive definition of disks 
that have holes. If we only consider objects with IR dips 
in their SEDs as transitional, the TD fraction is lower; 
including objects with decreasing emission at all wave- 
lengths increases the reported "tr ansitional disk" fractio n 
from 20% to 70% at ~10 Myr (IMuzerolle et al.l [2010h . 
This difference is not trivial and presents an unclear pic- 
ture when attempting to compare theoretical simulations 
of planet formation that predict disk holes with a "transi- 
tional disk" fraction that encompasses objects which do 
not have apparent evidence for cleared regions in their 
disks. 

Another related issue is encountered when trying to 
discern the disk clearing timescale. These timescales usu- 
ally include transitional and pre-transitional disks and 
evolved disks (which are optically thin), and are mea- 
sured with respect to full disks. Defining the boundaries 
between transitional, pre-transitional, evolved, and full 
disks is crucial in order to obtain an accurate estimate 
of the disk clearing timescale. TD are relatively easier to 
identify, whether looking at broad-band SEDs or colors. 
PTD are harder to tell apart from full disks based on col- 
ors alone. Evolved disks are also difficult to separate from 
full disks and this is the main driving force behind the dif- 
ferent disk clearing timescales reported in the literature. 



iCurrie fc Sicilia-Aguilarl ([20111 ) find a ~1 Myr timescale 
for inner disk clearing, which is longer than th e ~0.5 Myr 
timescale obtained by iLuhman et al.l ((2010). For the 
most part, both groups are using the same data. The 
main difference is what one considers a full disk versus 
evolve d (also called homologou sly depleted). As pointed 
out by iHernandez et al.l (|2010fl , the reported fraction of 
optically thin disks (i.e. evolved disks or homologously 
depleted d isks) is highly dependent on the cutoff. Obser- 
vationally, ICurrie fe Sicilia-Aguilarl ([201 If) use the lower 
quartile of Taurus. Therefore, by construction, 2 5% o f 
disks in Taurus are evolved disks. ILuhman et al.l ([20101 ) 
use a gap in the IR color-color diagrams which leads to 
fewer objects in this phase. 

Here we focus on the physical structure that could be 
underlying the observed SEDs by using physical models 
to motivate our interpretation. We find that it is possible 
to group objects based on their observed SEDs and find a 
general model-based interpretation to fit objects within 
a group. In our work we have cases of disks with holes 
and gaps and full disks (see Section |3~3|) . We emphasize, 
as noted previously in Section [4. 11 that possibly all disks 
have gaps that cannot be detected in SEDs. However, 
our goal here is to discern when one can identify a hole 
or gap in a disk based on its SED. 

In essence, all disks around TTS are "in transition" as 
they are all evolving in one way or another. The expec- 
tation is that all TTS with disks will eventually become 
diskless stars. However, they are not all necessarily going 
down the same evolutionary path. Here we suggest that 
the term "transitional" be used for a disk which appears 
to be undergoing a radical disturbance in the radial struc- 
ture of its inner disk (i.e., a hole or gap). While above we 
point out the inherent deficiencies in using the evolution- 
ary term "transitional" to define a disk, introducing new 
classification schemes to the literature is not warranted 
given the limitations of currently available data. 

One motivation for separating disks with holes and 
gaps from full disks is that it is not obvious these 
disks are undergoing the same type of disk clear- 
ing. Many researchers have suggested that the 
holes and gaps in disks obser ved to date are du e 
to planets (see discussion in Espai llat et al.l I2010T) . 
Simulations have shown that newly forming plan- 
ets will clear regions of the disk through ac c retion 
and tidal disturbance s (IGoldreich fe Tremaind 119801: 
Ward! 119881: iRice et all 120031: iPaardekooper fe Mellema j 
2004 HQuillen et al.ll2004J : iVarniere et al.ll2006l : iZhu et al.l 
1201 If ). It is less clear how a disk with weak emission at 
all wavelengths could be related to disk clea r ing ca used 
by a planet. As pointed out b y ICieza et ail (|2010f ) and 
ICurrie fe Sicilia-Aguilarl (|2011h . disks with weak MIR 
emission could instead be the result of another mech- 
anism that may have a different rate of evolution (e.g. 
photoevaporation). Alternatively, full disks with SEDs 
such as those in this paper could simply be the tail end 
of continuous distribution of full disks. This could be 
related to a large spread in disk properties (M, M<i, dust 
composition) in a given population as well as a distribu- 
tion in the initial conditions. Another possibility, that we 
cannot fully test in this paper due to a lack of millimeter 
detections, is that the full disks in our sample have expe- 
rienced a greater degree of dust settling than other full 
disks. In this case, we would see more of these disks in 
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older regions since settling is expected to increase with 
age. 

4.4. Mass Accretion Rates of Transitional and 
P re-transitional Disks 

Naii ta et al.1 (j2007j ) showed that the mass accretion 
rates of transitional disks in Taurus tend to be lower 
than those of full disks in the same region. If planets are 
the clearing agent in transitional disks, then lower mass 
accretion rates are expected onto the star since a giant 
planet that opens a gap in the disk will intercept and ac- 
crete material from the outer disk (jLubow k, D'Angelol 
2006) . To explore this further we compared the distribu- 
tion of mass accretion rates of full disks and transitional 
and pre-transitional disks in Taurus, Chamaeleon, and 
NGC 206 8 (Figure fl"2l s ee the Appendix for details). We 
note that iNajita et al.1 (|2007j ) used a broader definition 
of "transitional disk" which included all objects with less 
emission than the median SED of Taurus. As discussed 
in Section 14.31 the link between these disks and planet 
formation is less clear. Therefore, here we use our more 
restrictive definition of transitional and pre-transitional 
disks which includes only objects with holes and gaps. 
We find that the mass accretion rates of transitional and 
pre-transitional disks tend to be about 5 times lower than 
the full disks in these three regions. The median mass ac- 
cretion rate for the full disks isl.3xl0 -8 M© yr -1 and for 
the transitional and pre-transitional disks it is 3.1 xlO -9 
M Q yr" 1 . A Kolmogorov-Smirnov (KS) test indicates 
that the full disk and transitional and pre-transitional 
disk mass accretion rate samples are not drawn from the 
same distribution (the KS-probability is 0.02). 

While the mass accretion rates of transitional and pre- 
transitional disks are overall lower than those of full 
disks, they are still too high to be compatible with cur- 
rent models of disk clearing by planets. This is espe- 
cially seen in the cases of transitional and pre-transitional 
disks w ith higher mass a ccretion rates and large gaps and 
holes. iZhu et ail (|2011[) find that multiple planets are 
needed to open these large clearings in the dust distri- 
bution. However, more planets in the disk should lead 
to lower mass accretion rates onto the star than those 
observed. Our results suggest that possible planets in 
transitional and pre-transitional disks could be lowering 
the mass accretion rate onto the star somewhat, but that 
there is another mechanism taking effect that we have 
no t accounted for, possibly dust evolution as proposed 
bv IZhu eT al. (2011). More simulations of disk clearing 
by planets are needed to reconcile the large gap sizes and 
mass accretion rates currently observed. 

We also see that the transitional disks tend to have 
have lower mass accretion rates than the pre-transitional 
disks in our sample, by a factor of ~10. The median mass 
accretion rate for our five transitional disks is 9.7xl0~ 10 
Mq yr" 1 while the median for the ten pre-transitional 
disks in the sample is 8.8 xlO -9 M Q yr _1 . (More obser- 
vations of PTD and TD are needed to expand the sample 
size and confirm this result given that the KS-probability 
that the samples are drawn from the same distribution is 
0.27.) Given that the evolution and relationship between 
TD and PTD is not currently completely understood, 
the underlying reason for this apparent discrepancy in 
mass accretion rates is not obvious. One can speculate 
that the difference is due to the same mechanism clear- 



ing the holes and gaps in t h ese d isks. In the case of 
planet formation, IZhu et ahl ()2011[ ) find that the mass 
accretion rate onto the star will decrease with time as 
planets grow in the disk. The difference in mass accre- 
tion rate between PTD and TD could then possibly indi- 
cate that pre-transitional disks are in the early stages of 
planet formation while transitional disks are in the later 
stages. However, refinement of planet forming simula- 
tions is needed to study this further given the complex 
structure expected in the inner disk region. 

5. SUMMARY 

Here we modeled the broad-band SEDs of 15 disks in 
NGC 2068 and IC 348. We presented IRS spectra for 
all our targets as well as mass accretion rates estimated 
with U-band photometry obtained at the MDM Observa- 
tory and H a profiles from the MIKE spectrograph on the 
Magellan telescope. We also presented SMA millimeter 
data for some of our sources in IC 348. 

The observed SEDs of the objects in our sample are 
diverse, yet can be separated into three groups. Some of 
our targets have dips in both their NIR and MIR emis- 
sion, some have dips in only their MIR emission, and 
some have decreasing emission at all IRS wavelengths. 
We modeled the first group as transitional disks (i.e. ob- 
jects with holes in their disk's dust distribution), the sec- 
ond group as pre-transitional disks (i.e. objects with gaps 
in their disk's dust distribution), and the last group as 
full disks (i.e. objects with no cleared regions in their 
disks). We found that millimeter data are crucial in 
breaking model degeneracies between the amount of dust 
settling in the disk and the disk's mass. 

We discussed the limitations of the observations, 
namely that we currently do not have high enough reso- 
lution to discern very small gaps in disks, and the limita- 
tions of disk models, especially with respect to simulating 
the effects of dust settling and determining masses. We 
pointed out that much of the disagreement in the lit- 
erature over reported transitional disk frequencies and 
disk clearing timescales is mainly due to inconsistent 
application of the term "transitional" in the literature. 
We suggested that only objects showing evidence of an 
abrupt change in their radial disk structure be referred 
to as "transitional." Specifically, here we use "transi- 
tional disk" when referring to disks with holes and "pre- 
transitional disk" for disks with gaps. Finally, we com- 
pared the mass accretion rates of transitional and pre- 
transitional disks to full disks in Taurus, Chamaeleon, 
and NGC 2068 and find that PTD and TD have lower 
accretion rates overall. We also find that the TD have 
lower mass accretion rates than PTD, but due to our 
small sample more objects are needed to confirm this. 

Significant progress will be made in the near future on 
the issues raised in this paper. Herschel SPIRE will pro- 
vide us with a large, consistent sample of sub-millimeter 
fluxes to help break model degeneracies. With the high 
resolution of ALMA we can soon test if the above clas- 
sifications used in this paper hold and modify them if 
necessary. 
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referee for a careful review of the manuscript. C. E. was 
supported by the National Science Foundation under 



The Transitional Disk Class 



11 



Award No. 0901947. P. D. acknowledges a grant from 
PAPIIT-DGAPA UNAM. N. C. acknowledges support 
from NASA Origins Grant NNX08AFM 5154G. 

REFERENCES 

Andrews, S. M., & Williams, J. P. 2005, ApJ, 631, 1134 
Andrews, S. M., Wilner, D. J., Espaillat, C, Hughes, A. M., 

Dullemond, C. P., McClure, M. K., Qi, C, & Brown, J. M. 

2011, ApJ, 732, 42 
Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C, & 

Dullemond, C. P. 2009, ApJ, 700, 1502 
— . 2010, ApJ, 723, 1241 

Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651 

Barrado y Navascues, D., & Martin, E. L. 2003, AJ, 126, 2997 

Batalha, C. C, & Basri, G. 1993, ApJ, 412, 363 

Bernstein, R., Shectman, S. A., Gunnels, S. M., Mochnacki, S., & 
Athey, A. E. 2003, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Scries, Vol. 4841, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, ed. M. Iyc & A. F. M. Moorwood, 1694-1704 

Birnstiel, T., Ormel, C. W., & Dullemond, C. P. 2011, A&A, 525, 
A11+ 

Calvct, N., Briceno, C, Hernandez, J., Hoyer, S., Hartmann, L., 
Sicilia-Aguilar, A., Megeath, S. T., & DAlessio, P. 2005a, AJ, 
129, 935 

Calvet, N., D'Alessio, P., Hartmann, L., Wilner, D., Walsh, A., & 

Sitko, M. 2002, ApJ, 568, 1008 
Calvet, N., & Gullbring, E. 1998, ApJ, 509, 802 
Calvet, N., Patino, A., Magris, G. C, & DAlessio, P. 1991, ApJ, 

380, 617 

Calvet, N., et al. 2005b, ApJ, 630, L185 

Cieza, L. A., ct al. 2010, ApJ, 712, 925 

Cody, A. M., & Hillenbrand, L. A. 2010, ApJS, 191, 389 

Cohen, M., Megeath, S. T., Hammcrsley, P. L., Martm-Luis, F., 

& Stauffer, J. 2003, AJ, 125, 2645 
Currie, T., & Sicilia-Aguilar, A. 2011, ApJ, 732, 24 
D'Alessio, P., Calvet, N., & Hartmann, L. 2001, ApJ, 553, 321 
DAlessio, P., Calvet, N., Hartmann, L., Franco-Hernandez, R., & 

Servm, H. 2006, ApJ, 638, 314 
DAlessio, P., Calvet, N., Hartmann, L., Lizano, S., & Canto, J. 

1999, ApJ, 527, 893 
D'Alessio, P., Canto, J., Calvet, N., & Lizano, S. 1998, ApJ, 500, 

411 

D'Alessio, P., et al. 2005, ApJ, 621, 461 
Dodson-Robinson, S. E., & Salyk, C. 2011, ApJ, 738, 131 
Dullemond, C. P., & Dominik, C. 2004, A&A, 417, 159 
— . 2005, A&A, 434, 971 

Dullemond, C. P., Dominik, C, & Natta, A. 2001, ApJ, 560, 957 
Espaillat, C. 2009, PhD thesis, University of Michigan 
Espaillat, C, Calvet, N., D'Alessio, P., Hernandez, J., Qi, C, 

Hartmann, L., Furlan, E., & Watson, D. M. 2007a, ApJ, 670, 

L135 

Espaillat, C, Calvet, N., Luhman, K. L., Muzerolle, J., & 

D'Alessio, P. 2008a, ApJ, 682, L125 
Espaillat, C, Furlan, E., D'Alessio, P., Sargent, B., Nagcl, E., 

Calvet, N., Watson, D. M., & Muzerolle, J. 2011, ApJ, 728, 49 
Espaillat, C, et al. 2007b, ApJ, 664, Llll 
— . 2008b, ApJ, 689, L145 
— . 2010, ApJ, 717, 441 

Fischer, W., Edwards, S., Hillenbrand, L., & Kwan, J. 2011, ApJ, 
730, 73 

Flaherty, K. M., & Muzerolle, J. 2008, AJ, 135, 966 

Flaherty, K. M., Muzerolle, J., Rieke, G., Gutermuth, R., Balog, 

Z., Hcrbst, W., Megeath, S. T., & Kun, M. 2011, ApJ, 732, 83 
Flaherty, K. M., Muzerolle, J., Rieke, G., Gutermuth, R., Balog, 

Z., Herbst, W., Megeath, S. T., & Kun, M. 2012, ApJ, 

submitted 

Franchini, M., Morossi, C, & Malagnini, M. L. 1998, ApJ, 508, 
370 

Furlan, E., et al. 2006, ApJS, 165, 568 
— . 2009, ApJ, 703, 1964 
— . 2011, ApJS, 195, 3 

Goldreich, P., & Tremaine, S. 1980, ApJ, 241, 425 
Gullbring, E., Hartmann, L., Briceno, C, & Calvet, N. 1998, 
ApJ, 492, 323 



Haisch, Jr., K. E., Lada, E. A., & Lada, C. J. 2001, AJ, 121, 2065 
Hartigan, P., Hartmann, L., Kenyon, S., Hewett, R., & Stauffer, 

J. 1989, ApJS, 70, 899 
Hartmann, L., Calvet, N., Gullbring, E., & D'Alessio, P. 1998, 

ApJ, 495, 385 
Hcrbig, G. H. 1998, ApJ, 497, 736 

Hernandez, J., Morales-Calderon, M., Calvet, N., Hartmann, L., 
Muzerolle, J., Gutermuth, R., Luhman, K. L., & Stauffer, J. 

2010, ApJ, 722, 1226 

Hernandez, J., et al. 2007, ApJ, 671, 1784 
Higdon, S. J. U., et al. 2004, PASP, 116, 975 
Houck, J. R., et al. 2004, ApJS, 154, 18 

Houdebine, E. R., Mathioudakis, M., Doyle, J. G., & Foing, B. H. 

1996, A&A, 305, 209 
Huelamo, N., Lacour, S., Tuthill, P., et al. 2011, A&A, 528, L7 
Hughes, A. M., Wilner, D. J., Calvet, N., D'Alessio, P., Claussen, 

M. J., & Hogerheijde, M. R. 2007, ApJ, 664, 536 
Hughes, A. M., et al. 2009, ApJ, 698, 131 

Ingleby, L., Calvet, N., Hernandez, J., Briceno, C, Espaillat, C, 

Miller, J., Bergin, E., & Hartmann, L. 2011, AJ, 141, 127 
Ingleby, L., & et al. 2011, submitted to ApJ 

Isella, A., Carpenter, J. M., & Sargent, A. I. 2010, ApJ, 714, 1746 
Isella, A., & Natta, A. 2005, A&A, 438, 899 

Jensen, E. L. N., Mathieu, R. D., & Fuller, G. A. 1994, ApJ, 429, 
L29 

Jordi, K., Grebel, E. K., & Ammon, K. 2006, A&A, 460, 339 

Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 

Kraus, A. L., Ireland, M. J., Martinache, F., & Hillenbrand, L. A. 

2011, ApJ, 731, 8 

Kraus, A. L., & Ireland, M. J., ApJ, in press 
Lada, C. J., et al. 2006, AJ, 131, 1574 
Landolt, A. U. 1992, AJ, 104, 340 

Leboutciller, V., Bcrnard-Salas, J., Sloan, G. C, & Barry, D. J. 

2010, PASP, 122, 231 
Lima, G. H. R. A., Alencar, S. H. P., Calvet, N., Hartmann, L., & 

Muzerolle, J. 2010, A&A, 522, A104+ 
Lubow, S. H., & D'Angelo, G. 2006, ApJ, 641, 526 
Luhman, K. L., Allen, P. R., Espaillat, C, Hartmann, L., & 

Calvet, N. 2010, ApJS, 186, 111 
Luhman, K. L., Stauffer, J. R., Muench, A. A., Rieke, G. H., 

Lada, E. A., Bouvier, J., & Lada, C. J. 2003, ApJ, 593, 1093 
Mathis, J. S. 1990, ARA&A, 28, 37 

Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 

McClure, M. 2009, ApJ, 693, L81 

McClure, M. K., ct al. 2008, ApJ, 683, L187 

Merin, B., et al. 2010, ApJ, 718, 1200 

Miyake, K., & Nakagawa, Y. 1995, ApJ, 441, 361 

Muzerolle, J., Allen, L. E., Megeath, S. T., Hernandez, J., & 

Gutermuth, R. A. 2010, ApJ, 708, 1107 
Muzerolle, J., Calvet, N., & Hartmann, L. 1998, ApJ, 492, 743 
Muzerolle, J., Calvet, N., Hartmann, L., & D'Alessio, P. 2003, 

ApJ, 597, L149 

Najita, J. R., Strom, S. E., & Muzerolle, J. 2007, MNRAS, 378, 
369 

Natta, A., Testi, L., Muzerolle, J., Randich, S., Comeron, F., & 

Persi, P. 2004, A&A, 424, 603 
Olofsson, J., Benisty, M., Augcreau, J.-C, ct al. 2011, A&A, 528, 

L6 

Owen, J. E., Ercolano, B., & Clarke, C. J. 2011, MNRAS, 412, 13 
Paardekooper, S.-J., & Mellema, G. 2004, A&A, 425, L9 
Pott, J.-U., Perrin, M. D., Furlan, E., et al. 2010, ApJ, 710, 265 
Quillen, A. C, Blackman, E. G., Frank, A., & Varniere, P. 2004, 
ApJ, 612, L137 

Rice, W. K. M., Wood, K., Armitagc, P. J., Whitney, B. A., & 

Bjorkman, J. E. 2003, MNRAS, 342, 79 
Rigliaco, E., Natta, A., Randich, S., Testi, L., Covino, E., 

Herczeg, G., & Alcala, J. M. 2011, A&A, 526, L6+ 
Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 

2007, ApJS, 169, 328 
Sargent, B. A., Forrest, W. J., Tayrien, C, et al. 2009, ApJ, 690, 

1193 

Schaefer, G. H., Simon, M., Beck, T. L., Nelan, E., & Prato, L. 

2006, AJ, 132, 2618 
Sicilia-Aguilar, A., Henning, T., Dullemond, C. P., Patel, N., 

Juhasz, A., Bouwman, J., & Sturm, B. 2011, submitted to ApJ 
Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 



12 



Espaillat et al. 



Skrutskie, M. F., Dutkevitch, D., Strom, S. E., Edwards, S., 

Strom, K. M-, & Shure, M. A. 1990, AJ, 99, 1187 
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Strom, K. M., Strom, S. E., Edwards, S., Cabrit, S., & Skrutskie, 

M. F. 1989, AJ, 97, 1451 
Thalmann, C, et al. 2010, ApJ, 718, L87 
Uchida, Y., & Shibata, K. 1984, PASJ, 36, 105 
Varniere, P., Blackman, E. G., Frank, A., & Quillen, A. C. 2006, 

ApJ, 640, 1110 
Ward, W. R. 1988, Icarus, 73, 330 



Watson, D. M., et al. 2009, ApJS, 180, 84 

Weidenschilling, S. J., Spaute, D., Davis, D. R., Marzari, F., & 

Ohtsuki, K. 1997, Icarus, 128, 429 
Werner, M. W., et al. 2004, ApJS, 154, 1 
White, R. J., & Basri, G. 2003, ApJ, 582, 1109 
White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265 
White, R. J., & Hillenbrand, L. A. 2004, ApJ, 616, 998 
Yang, H., Johns-Krull, C. M., & Valenti, J. A. 2007, AJ, 133, 73 
Zhu, Z., Nelson, R. P., Hartmann, L., Espaillat, C, & Calvet, N. 

2011, ApJ, 729, 47 



APPENDIX 
DUST COMPOSITION OF SAMPLE 

We performed fitting of the silicate emission features visible in the IRS spectra and derived the mass fraction of 
amorphous and crystalline silicates in the disk (Tables |51 [7] and H]). See Ell for a discussion of the degeneracies inherent 
in deriving the dust composition. The results here should be taken as representative of a dust composition that can 
reasonably explain the observed silicate features in the SED. We leave it to future work to further constrain the mass 
fractions of silicates in these disks. 

For the inner wall of our pre-transitional objects, we adopted a silicate composition consisting solely of amorphous 
olivines. This is because the inner wall does not produce significant 10 /im silicate emission in most of the objects in 
this study and so we have no way to distinguish between pyroxene and olivine silicates in the inner wall. The exception 
is LRLL 31 where we need 60% amorphous olivine and 40% forsterite in the inner wall in order to fit the 10 /im silicate 
emission feature. For transitional and pre-transitional disks we changed the silicate composition in the optically thin 
dust region and outer wall to fit the SED (Tables [7J . For the full disks, we changed the silicate composition in the 
inner wall and disk (Table [6] and [8]) . The silicate composition was not allowed to vary between the inner wall and disk 
in the full disk models. 

COMMENTS ON FITTING ZZ TAU WITH A SETTLED, IRRADIATED ACCRETION DISK MODEL 

Here we present modeling of the SED of ZZ Tau using the disk model of D'Alcssi o et al.l (|2006l ) discussed in Sec- 
tion GUI Stellar parameters used in the disk model (T»=3470 K; L*=0.75 L Q ; M*=0.35 M Q ; R*=2.4 R Q ) were 
derived in the same manner as other objects in this work (see Section [OJ us ing a spectral type of M3 adopted from 
iKenvon fc Hartma nn (1995) and a visual extinction (Ay) of 0.98 taken from iFurlan et ail ((2006). We adopt a mass 
accretion rate of 1.3xl0~ 9 M Q yr" 1 from lWhite fc Ghezl (|2001lh We note that this is slightly higher than the mass 
accretion rate measured from U-band photometry (9xl0 -10 M Q yr _1 ). 

In Figure [T2] we present two models with differe nt outer radii. In one model we use an outer disk radius of 100 AU for 
comparison with previous modeling performed bv lCurrie fc Sicilia-Aguilarl f|201 ID . The best-fit parameters are e=0.001 
and a=0.02 and this disk has a mass of 5xl0 -4 M . We also explored disks with other e values. A disk with a higher 
e of .01 needs an a of 0.2 t o the fit the SED, but this a results in a viscous timescale shorter than the lifetime of the 
disk partmann et al.lll998D . A disk with e=0.0001 and a=0.002 with a higher mass of 0.005 Mq is excluded by the 
millimeter upper limits. Here we use amorphous silica tes to fit the disk with a ma;c =10 /zm. Crystalline silicate features 
are evident in the IRS spectrum (|Sargent et al.l 120091) . but we leave a detailed fit to McClure et al. (in preparation). 
In Figure [T4l we show that ZZ Tau is optically thick to its own radiation (tr oss >1) o ut to about ~1 AU in the disk. 
Over 80% of the emission seen at 40 /jm is from within these radii (see Figure 2.16 in lEspaillati[2009l ). Therefore, the 
optically thick part of the disk of ZZ Tau dominates the emission seen in the IRS spectrum. 

We also present a model with an outer radius of 3 AU in Figure [T31 This is because ZZ Tau is a close binary with a 
separation of 0.06" ([Schaefer et al.l 12006). which corresponds to 8 AU at the distance of Taurus (140 pc). Therefore, 
the IRS spectrum presented here includes both o bjects. (ZZ Tau IRS, which is 36"away, did not enter the IRS slit and 
could be a wide companion (|Furlan et al.l [20TlT ) If a circumbinary disk is present, its inner edge would be located 
at ^16 AU according to expectations of dynamical clearing by companions (jArtvmowicz fc LubowllT994[) . However, 
we detect NIR blackbody emission which indicates tha t instead we are seeing a cir cumprimary disk. In this case, the 
outer edge of the disk would be truncated to ^3 AU (jArtymowicz fc Lubowl P~994[ ) and have a mass of 2xl0 -5 M . 
As discussed in Section lB. 3.21 and shown in Figure Q3J since the IRS emission is dominated by the inner AU of the disk, 
changing the outer radius of the disk does not significantly alter the IR emission. 

COMMENTS ON THE DISTRIBUTION OF MASS ACCRETION RATES 

When plotting the distribution of mass accretion rates of full disks and transitional and pre-transitional disks in 
Taurus, Chamaeleon, and NGC 2068 (Figure IT2")) we restricted ourselves to transitional and pre-transitional disks 
whose SEDs have be en modeled. Mass accretion rates for transitional and pre-transitional disks are taken from 
lEspaillat et al.l (|2011[) for Tauru s and Chamaeleon and from this work for NGC 2068. The mass accretion rates for full 
disks in Taurus were taken from Naii ta et al.l (|2007f) . Since here we use a different definition of "transitional disk" than 
used in that work, we use the mass accretion rates for objects that do not overlap wit h what we label a pre-tra nsitional 
or transitional disk. For Chamaeleon, mass accretio n rates for full disks are fro m lHartmann et al.l ()1998l ) and for 
NGC 2068 mass accretion rates for full disks are from iFlahertv fc Muzerollel (|2008l ). We do not include objects which 
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have upper limits for their mass accretion rates or those that are known to be binaries. We also excluded IC 348 in 
this analysis since, to the best of our knowledge, there are no mass accretion rates in the literature for the full disks 
in this region. In addition, IC 348 is older than Taurus, Cham aeleon, and NGC 20 68 which may bias the results given 
that mass accretion rates are known to decrease with age (e.g. iCalvet et alJl2005al ). In total we have 45 full disks and 
15 transitional and pre-transitional disks. 

We note that the maj ority of these mass accretion rates are derived using U-band photometry and the relation 
in iGul lbring et al. (1998J). The exceptions are objects in NGC 2068. The mass accretion rates for transitional and 
pre-tr ansitional disks in NGC 206 8 are taken from this work and the mass accretion rates for the full disks are adopted 
from (jFlahertv fc Muzerolldl2008l ). In Section [3XTJ we discuss the derivation methods used in both works. In short, 
the typical error (a factor of 2-3) inherent to mass accretion rate estimation methods should not lead to systematic 
differences between different samples. 
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Table 1 

Target Sample 



Name 


Region 


RA 


DEC 


FM 177 


NGC 2068 


05h45m42s 


-00dl2m05s 


FM 281 


NGC 2068 


05h45m53s 


-00dl3m25s 


FM 515 


NGC 2068 


05h46ml2s 


+00d32m26s 


r iVi ool 


INO^ zUOo 


nt^Vi A ft m 1 On 

Uon4onilys 


-UUqUOIIIooS 


FM 618 


NGC 2068 


05h46m23s 


-00d08m53s 


LRLL 2 


IC 348 


03h44m35s 


+32dl0m04s 


LRLL 6 


IC 348 


03h44m37s 


+32d06m45s 


LRLL 21 


IC 348 


03h44m56s 


+32d09ml5s 


LRLL 31 


IC 348 


03h44ml8s 


+32d04m57s 


LRLL 37 


IC 348 


03h44m38s 


+32d03m29s 


LRLL 55 


IC 348 


03h44m31s 


+32d00ml4s 


LRLL 67 


IC 348 


03h43m45s 


+32d08ml7s 


LRLL 68 


IC 348 


03h44m29s 


+31d59m54s 


LRLL 72 


IC 348 


03h44m23s 


+32d01m53s 


LRLL 133 


IC 348 


03h44m42s 


+32dl2m02s 



Note. — Target IP's are taken from 

IFlahertv fc Muzeroflg <2008D and ILuhman et al.l d2003ft 
for targets in NGC 2068 and IC 348, respectively. 



Table 2 

IC 348 Optical Photometry 



Target 


V 


U-V 


V-R 


V-I 


LRLL 


2 


sat. 


sat. 


sat. 


sat. 


LRLL 


6 


sat. 


sat. 


sat. 


sat. 


LRLL 


21 


15.68±0.03 


3.51±0.06 


1.35±0.08 


sat. 


LRLL 


31 


19.30±0.08 


<0.8 


2.00±0.09 


3.89±0.09 


LRLL 


37 


15.85±0.04 


2.68±0.25 


1.24±0.05 


1.19±0.06 


LRLL 


55 


21.68±0.33 


<-1.58 


2.03±0.22 


4.19±0.38 


LRLL 


67 


16.23±0.02 


2.04±0.11 


1.16±0.04 


2.50±0.04 


LRLL 


68 


17.49±0.03 


2.53±0.31 


1.64±0.03 


3.46±0.03 


LRLL 


72 


17.55±0.03 


2.38±0.20 


1.62±0.03 


3.32±0.03 


LRLL 


133 


20.00±0.30 


<0.1 


1.50±0.50 


3.70±0.30 



Note. — We use "sat." to refer to observations that were sat- 
urated and note upper limits for bands in which sources were not 
detected. 



Table 3 

Source Properties 



Target 


Ay 


Spectral 


T* 


L» 


M* 




Al 


Lace 


M 






Type 


(K) 


(M ) 


(M Q ) 


(Ro) 


( 10- 8 M yr" 1 ) 


(Lo) 


Source 


FM 177 


1.6 


K4 


4590 


1.0 


1.2 


1.5 


0.004 


0.0009 


Ha 


FM 281 


2.0 


Ml 


3720 


0.4 


0.5 


1.6 


0.002 


0.0002 


Ha 


FM 515 


1.6 


K2 


4900 


2.5 


1.5 


2.2 


3.10 


0.68 


Ha 


FM 581 


4.1 


K4 


4590 


4.1 


1.6 


3.1 


2.57 


0.40 


Ha 


FM 618 


2.9 


Kl 


5080 


2.2 


1.5 


1.9 


1.21 


0.29 


Ha 


LRLL 2 


3.8 


A2 


8970 


57.1 


2.8 


3.1 








LRLL 6 


3.9 


G3 


5830 


16.6 


2.4 


4.0 








LRLL 21 


4.7 


K0 


5250 


3.8 


1.6 


2.4 


0.20 


0.04 


Ha 


LRLL 31 


8.6 


G6 


5700 


5.0 


1.6 


2.3 


1.4 


0.3 


Fll 


LRLL 37 


2.8 


K6 


4205 


1.3 


0.9 


2.2 


0.13 


0.02 


U-band 


LRLL 55 


8.5 


M0. 5 


3850 


1.0 


0.6 


2.2 








LRLL 67 


2.0 


M0. 75 


3720 


0.5 


0.5 


1.8 


0.01 


0.001 


Ha 


LRLL 68 


2.1 


M3.5 


3470 


0.5 


0.3 


2.0 


0.04 


0.002 


U-band 


LRLL 72 


3.0 


M2.5 


3580 


0.7 


0.1 


2.1 


<0.0003 


<0.00001 


Ha 


LRLL 133 


3.6 


M5 


3240 


0.2 


0.2 


1.5 


<0.8 


<0.78 


U-band 



Note. — Spectral types for objects in NGC 2068 and IC 348 arc adopted from Flahcrtv & Muzcrollc 
(2008) and iLu hrnanet ahM2 003l) , respectively, except in the case of LRLL 31 where we adopt the spectral 
type of Flaherty ct al. (2011). T* is taken from Kcnyon Sz Hartmann (1995), based on the adopted spectral 
type. L„, M* and R* are calculated in this work. measuremen ts for most of the obj ects arc from this 

work except for LRLL 2 and LRLL 6 where this value is adopted from Luhman ct al. (2003). The last column 
lists the me thod with which our mass accretion rates were calculated: Ha and U-band are from this work; 
Fll is from Flaherty ct al. (2011). For sources where mass accretion rate estimates arc not listed here, we 
adopt lxlO -8 M yr -1 . 
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Table 4 

Wall Properties of Sample 



_L OjI >^*" u 


Disk 






Wall 






Outer Wall 




(!) 


lype 
(2) 


&max 

(c m ) 

(3) 


wall 

( K ) 

(A \ 

(4) 


(AU) 
(5) 


Dl 

(AU) 
(6) 


&max 

(pm) 
(7) 


wall 

(K) 
(8) 


/ tea 1 1 

(AU) 
(9) 


DO 

(AU) 

/-in\ 

(10) 


b M 177 


rnl S 

1 D 










0.25 


90 


3.1 


49 


FM 281 


TD 


_ 




~ 




0.25 


90 


.3.6 


31 


LRLL 67 


TD 










5.0 




1.0 


10 


LRLL 72 


TD 










1.0 


190 


0.6 


5 


LRLL 133 


TD 










0.25 


180 


0.6 


i 


FM 515 


PTD 


10 


1700 a 


0.0071 


0.12 


0.25 


150 


5 


45 


FM 618 


PTD 


1.0 


1400 


0.0055 


0.22 


5.0 


180 


0.7 


11 


LRLL 21 


PTD 


2 


1800 6 


0.0017 


0.13 


2.0 


220 


0.9 


9 


LRLL 31 


PTD 


1.0 


1400 


0.01 


0.32 


5.0 


180 


1.5 


11 


LRLL 37 


PTD 


0.25 


1400 


0.0098 


0.17 


0.25 


240 


0.6 


5 


FM 581 


FD 


0.25 


1400 


0.016 


0.3 










LRLL 2 


FD 


0.25 


1400 


0.013 


1.68 










LRLL 6 


FD 


1.0 


1400 


0.005 


0.54 










LRLL 55 


FD 


0.25 


1400 


0.02 


0.14 










LRLL 68 


FD 


10.0 


1400 


0.0023 


0.07 











Note. — Col (1): Name of target. Col (2): Assigned classification for our targets. Wc label objects as 
transitional disks (TD), prc-transitional disks (PTD), and full disks (FD). Col (3): Maximum grain size 
of dust used for the inner wall of the disk. The superscript i denotes "inner wall." Col (4): Temperature 
of the inner wall. Col (5): Height of the inner wall. Col (6): Radius of the inner wall. Col (7): Maximum 
grain size of dust used for the outer wall of the disk. The superscript o denotes "outer wall." Col (8): 
Temperature of the outer wall. Col (9): Height of the outer wall. Col (10): Radius of the outer wall. 
a For FM 515, the inner wall model required a temperature of 1700 K in order to fit the slope of the 
NIR emission. In some cases, temperatures >1400 K for the inner wall have also been needed to fit the 
SED previously (e.g. T35, UX Tau A; lEspaillat eTaI|[20To| , Ell). 

k We adopt a temperature of 1800 K for the inner wall based upon NIR SpeX spectral fitting 
IFlaherty eT"aH J2012T ). 



Table 5 

Model Runs for LRLL 68 



Model 


% 


Hd 

(AU) 


e 


a 


M-disk 

(M«) 


1 


60 


100 


0.0001 


0.0006 


0.005 


2 


60 


100 


0.001 


0.006 


0.0005 


3 


60 


100 


0.01 


0.06 


0.00004 


4 


60 


300 


0.001 


0.006 


0.001 


5 


60 


20 


0.001 


0.006 


0.0001 


6 


20 


100 


0.001 


0.006 


0.001 


7 


40 


100 


0.001 


0.006 


0.001 


8 


80 


100 


0.001 


0.006 


0.001 



Table 6 

Mass Fraction (in %) of Silicates in Outer Wall in Transitional and 
Pre-transitional Objects 



Target 


Amorphous 


Amorphous 


Crystalline 


Crystalline 


Crystalline 




Olivine 


Pyroxene 


Forsterite 


Enstatite 


Silica 


FM 177 


80 




10 




10 


FM 281 


90 








10 


FM 515 


100 










FM 618 


100 










LRLL 21 


80 




20 






LRLL 31 


80 






20 




LRLL 37 


70 




20 


10 




LRLL 67 


100 










LRLL 72 


60 


30 


5 


5 




LRLL 133 


90 




10 
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Table 7 

Mass Fraction (in %) of Silicates in Full Disks 



Target 


Amorphous 


Amorphous 


Crystalline 


Crystalline 


Crystalline 




Olivine 


Pyroxene 


Forstcritc 


Enstatitc 


Silica 


FM 581 


60 




20 


10 


10 


LRLL 2 


60 




20 


10 


10 


LRLL 6 


50 




20 


10 


20 


LRLL 55 


80 




10 


5 


5 


LRLL 68 


85 




5 


5 


5 



Table 8 

Mass Fraction (in %) of Silicates in Optically Thin Dust Region in Transitional 
and Pre-transitional Objects 



Target 


Amorphous 


Amorphous 


Crystalline 


Crystalline Crystalline 




Olivine 


Pyroxene 


Forstcritc 


Enstatitc 


Silica 


FM 177 


20 




40 


40 




FM 281 


80 




20 






FM 515 


100 










FM 618 


100 










LRLL 21 


50 




20 




30 


LRLL 31 












LRLL 37 


40 




20 


20 


20 


LRLL 67 


95 




5 






LRLL 72 












LRLL 133 


90 




10 






Note. - 


For LRLL 31, 


we did not need 


an optically thin region 


to account for 



the 10 fim silicate emission. This feature comes from the inner wall which has 60% 
amorphous olivine and 40% crystalline forstcrite. 



Table 9 

Properties of Optically Thin Dust Region in Transitional and Pre-transitional Objects 



Target 


Organics (%) 


Troilitc (%) 


Silicates (%) 


M dust (10-^M Q ) 




FM 177 


42 


11 


47 


2 


5 


FM 281 


19 


15 


66 


0.9 


0.25 


FM 515 


19 


15 


66 


0.8 


0.25 


FM 618 


12 


9 


79 


2 


3 


LRLL 21 


19 


15 


66 


0.6 


2 


LRLL 31 












LRLL 37 


26 


14 


60 


7 


5 


LRLL 67 


42 


11 


47 


6 


5 


LRLL 72 












LRLL 133 


42 


11 


47 


6 


0.25 



Note. - 



We use use a minimum grain size of 0.005 fim in the optically thin dust region. 
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400 



Figure 1. High-resolution H a line profiles obtained using the MIKE spectrograph for three of our IC 348 targets and all five of our 
NGC 2068 targets. Wide and asymmetric profiles (i.e. FM 515, FM 581, FM 618, LRLL 21) in dicate substantial acc retion onto the star 
and gas in the inner disk while narrow profiles indicate low accretion rates onto the star (e.g. White & Basri 2003). We note that the 
spectrum of LRLL 72 has been binned up by a factor of ~5 here for clarity. 
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Figure 2. SEDs of transitional disks in our sample. We show both the observed fluxes (red) and dereddened fluxes (bl u e; see Table [3] for 
Ay). In the NIR the emission is similar to that of the stellar photosphere (broken magenta line; Kcnvon & Hartmann (1995)), but rises 
in the MIR and at longer wavelengths. This indicates that the dust in the inner disk has been removed and that there is a hole in the 
disk. Open circles correspond to ground-based U-, B-, V-, R-, I-, and L-b and photometry; closed circles are 2MASS J-, H-, and K-band 
photometry; triangles are Spitzer IRAC and MIPS photometry. See Section l3.1l for data references. Solid lines are Spitzer IRS spectra from 
this work. (A color version of this figure is available in the online journal.) 
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Figure 3. SEDs of pre-transitional disks in our sample. For LRLL 31, FM 515, and FM 618, there is significant NIR emission above 
the stellar photosphere, but there is also a dip in the MIR emission and the flux increases at longer wavelengths. This indicates a gap at 
intermediate radii in the disk. In the case of LRLL 37, there is no clear "dip" in the emission; this object was cl assif ied as a pre-transitional 
disk based on its strong silicate emission which could not be re produced with a full disk model (see Section I3.3D ■ We classify LRLL 21 
as a PTD here even though its NIR excess is weak in this epoch; Flaherty ct al. (2012) find that LRLL 21 has significant NIR blackbody 
emission in more recent iRS observations. The color scheme, symbols, and lines are the same as in Figure[2] (A color version of this figure 
is available in the online journal.) 
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Figure 4. SEDs of full disks in our sample. The emission decreases at all wavelengths. The color scheme, symbols, and lines are the 
same as those used in Figure [2] (A color version of this figure is available in the online journal.) 
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Figure 5. IR SEDs and disk models_(solid, black lines) of tr ansitional and pre-transitional disks in our sample. Here we show only the 
dereddened data from Figures [2] and [3] Refer to Section 13.31 and Table [4] for model details. Separate model components are the stellar 
photosphere (magenta dotted line), the inner wall (gray short-long-dashed line), the outer wall (red dot-short-dashed line), and the optically 
thin small dust located within the inner disk (green long-dashed line). (A color version of this figure is available in the online journal.) 
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Figure 6. Left: IR SED and full disk model fit of LRLL 37. We cannot successfully reproduce the 10 /^m and 20 /an emission. This 
disk model has e=0.001, ct=0.0008, and a m ax=0.25fira. Right: IR SED and pre-transitional disk model of LRLL 37. We can fit the strong 
silicate emission with a model that has a gap in the disk with most of the silicate emission arising from small, optically thin dust within 
the gap. The color scheme, symbols, and lines are the same as in Figure[2] (A color version of this figure is available in the online journal.) 
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Figure 7. Left: IR SEDs and disk models of the pre-transitional disk LRLL 31 and the transitional disk LRLL 67. Right: Broad-band 
SEDs and disk models of LRLL 31 and LRLL 67 which include SMA millimeter data (from this work). With millimeter data we can 
constrain the outer disk and so we include both the outer wall and the disk behind it in the model presented here (red dot-short-dash). 
We also show the contribution to the SED from the inner wall of LRLL 31 (gray short-long-dash) and the optically thin small dust located 
within the inner hole of LRLL 67 (green long-dash). Other symbols and lines are the same as used in Figure [2] (A color version of this 
figure is available in the online journal.) 
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Figure 8. IR SEDs and disk models of full disks in our sample. We can explain all of these disks with irradiated accretion disk models 
which incorporate dust settli ng. I n the case of FM 581, we need to truncate the outer disk radius to <1 AU in order to reproduce its very 
steep slope. Refer to Section 13.31 and Table [4] for model details. Here we show the contribution to the SED from the stellar photosphere 
(magenta dotted line), the inner wall (gray short-long-dash), and the outer disk (red dot-short-dash). (A color version of this figure is 
available in the online journal.) 
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Figure 9. Disk model fits to the SED of LRLL 68 with e/o=6. Models with the same e-to-o ratio will have similar emission in the IR 
but substantially different emission in the millimeter. Here we show the following models from Table[5] Model 1 (e=0.0001; short-dashed 
line), Model 2 (e=0.001; dotted line), and Model 3 (e=0.01; long-dashed line). Note that each model falls below the millimeter upper limit 
presented in this work (open triangle). This reflects the need for millimeter detections to constrain disk models. (A color version of this 
figure is available in the online journal.) 
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1 10 100 1000 

X(fim) 

Figure 10. Disk models of LRLL 68 with different R<j, keeping e and a fixed. Varyingthe outer radius of the disk will lead to changes in 
the disk mass, but small changes in the millimeter flux. Models shown are from Table \S\ Model 4 (R d =300 AU; short- long-dashed line), 
Model 2 (R£j=100 AU; dotted line), and Model 5 (R ( ;=20 AU; dot- long-dashed line). (A color version of this figure is available in the 
online journal.) 
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Figure 11. Disk models of LRLL 68 with different i. We find that changing the inclination of the disk, while keeping all other parameters 
fixed, will lead to substantial differences in the IR emission, but will not significantly affect the millimeter flux. In order of increasing 
inclination (i.e. from nearly face-on to nearly edge-on) we show the following models from Table [5] Model 6 (i=20°; short-long-dashcd 
line), Model 7 (i=40°; dot- long-dashed line), Model 4 (i=60°; dotted line), Model 8 (i=80°; dot- short-dashed line). (A color version of 
this figure is available in the online journal.) 
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Figure 12. Distribution of mass accretion rates for Taurus, Chamaeleon, and NGC 2068. We separate full disks (white area) and 
transitional and pre-transitional disks (shaded area); overall, transitional and pre-transitional disks have lower mass accretion rates than 
full disks. (A color version of this figure is available in the online journal.) 
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Figure 13. Models of ZZ Tau using an irradiated accretion disk with dust settling. We show a disk with amorphous silicates and an 
outer radius of 3 AU (solid) an d 100 AU (broken). The col or scheme, symbols, and lines are t he same as in F i gur es [2] and [8] U-,B-, V-, 
R-, and I- photomet ry are fr om Kcnyon & Hartmann (1995j); J-, H-, and K-band data are fro m Skrutskic ct al. (2006); Spitzer IRAC and 
MIPS data are from lLuhman et alJ 112 0101) and the I RS spectrum is from Furlan et al. (2006). All millimeter upper limits are taken from 
Jensen ct al. ( 1994) and Andrews & Williams (2005). (A color version of this figure is available in the online journal.) 
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Figure 14. Rosseland mean optical depth of ZZ Tau for the model shown in Figure [13] ZZ Tau is optically thick to its own radiation 
(tHoss>1) out to ~1 AU in the disk. It is these innermost disk radii that dominate the emission seen in the IRS spectrum. 



